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ABSTRACT 

We present a detailed analysis of Swift multi- wavelength observations of GRB 070110 and its re- 
markable afterglow. The early X-ray light curve, interpreted as the tail of the prompt emission, 
displays a spectral evolution already seen in other gamma-ray bursts. The optical afterglow shows a 
shallow decay up to ^2 d after the burst, which is not consistent with standard afterglow models. The 
most intriguing feature is a very steep decay in the X-ray flux at ~2xl0 s after the burst, ending an 
apparent plateau. The abrupt drop of the X-ray light curve rules out an external shock as the origin 
of the plateau in this burst and implies long-lasting activity of the central engine. The temporal and 
spectral properties of the plateau phase point towards a continuous central engine emission rather 
than the episodic emission of X-ray flares. We suggest that the observed X-ray plateau is powered by 
a spinning down central engine, possibly a millisecond pulsar, which dissipates energy at an internal 
radius before depositing energy into the external shock. 
Subject headings: gamma rays: bursts; X-rays: individual (GRB 070110) 



1. INTRODUCTION 

The Swift Gamma Ray Burst Explorer (Gehrels et al. 
2004) is a multi-wavelength observatory specifically de- 
signed to study gamma-ray burst (GRB) evolution from 
their early stages. It is equipped with a wide-field instru- 
ment, the Burst Alert Telescope (BAT; Barthelmy et al. 
2005a), covering the 15-350 keV energy band, and two 
narrow field instruments, the X-Ray Telescope (XRT; 
Burrows et al. 2005a) and the Ultraviolet/Optical Tele- 
scope (UVOT; Roming et al. 2005), covering the 0.2-10 
keV band and the 1700-6500 A wavelength range, respec- 
tively. 

During the first two years of its mission Swift detected 
^200 bursts, providing well-sampled X-ray and optical 
afterglow light curves. Nousek et al. (2006) and Zhang 
et al. (2006a) identified a common pattern in the X-ray 

1 INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica, 
Sezione di Palermo, via Ugo la Malfa 153, 90146 Palermo, Italy 

2 Department of Physics and Astronomy, University of Leicester, 
Leicester, LEI 7RH, UK 

3 Dipartimento di Scienze Fisiche ed Astronomiche, Sezione di 
Astronomia, Universita di Palermo, Piazza del Parlamento 1, 90134 
Palermo, Italy 

4 Department of Physics and Astronomy, University of Nevada 
Las Vegas, NV, USA 

5 Universita degli studi di Milano-Bicocca, Dipartimento di 
Fisica, piazza delle Scienze 3, 1-20126 Milano, Italy 

6 INAF - Osservatorio Astronomico di Brera, via Emilio Bianchi 
46, 1-23807 Merate (LC), Italy 

7 NASA, Goddard Space Flight Center, Greenbelt, MD 20771, 
USA 

8 Department of Astronomy & Astrophysics, 525 Davey Lab., 
Pennsylvania State University, University Park, PA 16802, USA 

9 Universities Space Research Association, 10211 Wincopin Cir- 
cle, Suite 500, Columbia, MD 20144 USA 

10 Department of Physics, Guangxi University, Nanning 530004, 
China 

11 National Astronomical Observatories/ Yunnan Observatory, 
CAS, Kunming 650011, China 

12 UCL Mullard Space Science Laboratory, Holmbury St Mary, 
Dorking, Surrey RH5 6NT 



light curves, described as consisting of three power law 
segments: an early steep decay with a temporal slope 
3<a<5 lasting up to ^300 s, followed by a shallower 
phase (0.5<a<l); the slope of the light curve steepens 
again (1.0<a<1.5) at ~10 3 -10 4 s after the burst onset. 
X-ray flares, overlaid on the underlying power law decay, 
have been detected up to 10 5 s after the BAT trigger 
(Burrows et al. 2005b; Chincarini et al. 2007). 

O'Brien et al. (2006) modeled the observed shape of 
early X-ray light curves with an exponential decay that 
relaxes into a power law and, in most cases, presents 
the shallow phase. They showed that the initial rapid 
decay is a smooth extension of the prompt emission, 
probably due to emission from large angles relative to 
the observer's line of sight (Kumar & Panaitescu 2000). 
A long-lasting energy injection into the forward shock, 
due either to a late internal activity or to a radial dis- 
tribution of Lorentz factors, has been invoked to inter- 
pret the shallower stage of the X-ray light curve (Rees & 
Meszaros 1998; Zhang & Meszaros 2002b; Granot & Ku- 
mar 2006). This explanation is consistent with the final 
smooth steepening into the "standard afterglow" decay 
(Meszaros & Rees 1997). 

The afterglow of GRB 070110 cannot be easily traced 
back to this well-known scenario. Its early X-ray light 
curve seems to display a canonical shape, with an initial 
steep decay and then a very flat plateau phase (a ~0.05), 
but ^2xl0 4 s after the trigger the count rate drops 
abruptly by more than one order of magnitude with a 
slope a>7 (Sbarufatti et al. 2007; Krimm et al. 2007a). 
The behavior of the optical afterglow seen by the UVOT 
is very different, showing a shallow smooth decay remain- 
ing fairly bright at later times. 

In the light of the unique properties of its afterglow, the 
Swift team identified GRB 070110 as a 'burst of interest' 
encouraging a follow-up campaign (Krimm et al. 2007b). 

The optical afterglow of GRB 070110 (Krimm et al. 
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2007a) was also detected by the ESO VLT equipped with 
the FORS2 instrument. A redshift of z=2.352±0.001 has 
been inferred on the basis of several absorption features 
in the spectra (Jaunsen et al. 2007). Further VLT obser- 
vations, performed ~10.7 d after the burst, detected a 
fainter than predicted afterglow, suggesting the presence 
of an optical break (Malesani et al. 2007a). 

In this paper we report on the 7-ray, X-ray and optical 
observations performed by Swift. The paper is organized 
as follows: in § 2 we present a multiwavelength timing 
and spectral analysis of both the prompt and the after- 
glow emission; in § 3 we discuss our results. Finally, in 
§ 4 we summarize our findings and conclusions. 

Throughout the paper, times are given relative to 
the BAT trigger time To, t=T-To, and the convention 
F„ jt °c v~P\r a has been followed, where the energy in- 
dex (3 is related to the photon index T = (3 + 1. We have 
adopted a standard cosmology model with Hubble con- 
stant Ho=70kms~ 1 Mpc -1 and cosmological constants 
n A = 0.73, fi M = 0.27 (Spergel et al. 2006). The phe- 
nomenology of the burst is presented in the observer time 
frame, unless otherwise stated. 

All the quoted errors are given at 90% confidence 
level for one interesting parameter (Ax 2 =2.706, Lamp- 
ton et al. 1976). 

2. DATA ANALYSIS 
2.1. Observations 

GRB 070110 triggered the Swift BAT at 07:22:41 UT 
on 10th January, 2007. The Swift narrow field in- 
struments, XRT and UVOT, began observing 93 s and 
85 s after the trigger, respectively. An accurate after- 
glow position was rapidly determined by the UVOT at 
R.A.=00 h 03 m 39.20 s , Dec.=-52°58'26.3" (J2000, Krimm 
et al. 2007a), with an uncertainty of 1 arcsec. 

XRT observations began with an initial 2.5 s im- 
age mode frame and then, as the source was bright 
(~40 ctss -1 ), collecting data in Window Timing (WT) 
mode. The XRT automatically switched to Photon 
Counting (PC) mode when the source decreased to 
~2 ctss -1 . Follow-up observations lasted 26 days for 
a total net exposure of 165 s in WT mode and 330 ks in 
PC mode. 

The UVOT took a short exposure with the V filter 
while the spacecraft was settling at the end of the initial 
slew. This exposure was followed by a "finding chart" 
exposure with the White filter lasting 100 s and then an 
exposure with the V filter lasting 400 s. UVOT then 
began its usual procedure of cycling through its 3 visible 
filters {V, B, and U) and 3 UV filters (UVW1, UVM2, 
and UVW2). The optical afterglow was detected in the 
White, U, B, and V filters, but not in the UV filters. 
The lack of detection of the afterglow in the UV filters 
is consistent with the measured redshift. A total of 219 
UVOT exposures were taken in the first 6.7 days, after 
which the afterglow fell below the detection thresholds. 

2.2. Gamma-ray data 

We analyzed BAT event data using the standard BAT 
analysis software included in the NASA's HEASARC 
software (HEASOFT, version 6.1.2). Fig. 1 presents the 
BAT mask-weighted light curve in the 15-150 keV en- 
ergy band. It shows a first main peak at t~0 s and then 
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Fig. 1.— BAT Mask-weighted light curve (15-150 keV) of 
GRB 070110. Dashed vertical lines mark the T90 duration. 

a decay on which several peaks are superposed; emission 
is visible until t~100 s. We estimated the burst dura- 
tion, defined as the interval containing 90% of the total 
observed fluence, to be T 90 (15-150 keV)=89±7 s. 

Spectra were created using the task batbinevt, updat- 
ing relevant keywords with batupdatephakw. The corre- 
sponding response matrices were generated by the task 
batdrmgen. Systematic errors were properly added to 
the spectra. Since the spacecraft slew started ^40 s af- 
ter the trigger, we created two different spectra (before 
and during the slew) with the appropriate response ma- 
trices and performed a joint fit. 

The time-averaged T90 spectrum (from -0.4 s to 89 s) 
can be fitted with a simple power law of photon in- 
dex 1.57±0.12. A cut-off power law or a Band model 
(Band et al. 1993) do not provide a better description 
and cannot constrain the peak energy value. In the 
hardness ratio light curves, comparing different BAT en- 
ergy bands, no sign of a significant spectral evolution is 
present throughout the prompt emission. 

The fluence over the 15-150 keV band is 
(L.8^q' 3 )x10 -6 erg cm -2 , from which we derive an 
observed isotropic energy of 2.3xl0 52 erg. This value 
can be considered as a lower limit to the isotropic 
energy, E 7j ; so , which is defined over a larger energy band 
(1 keV-10 McV in the source rest frame). Extrapolating 
our best fit model over the whole 1 keV-10 MeV 
(source rest frame) we estimate an upper limit of 
E 7) i SO <1.3x 10 53 erg. This value is derived under the 
extreme assumption that the best description of the 
spectrum is Band law (Band et al. 1993) with a peak 
energy E p >10 MeV in the source rest frame. 

2.3. X-ray data 

XRT data were processed using the XRTDAS software 
package (v. 2.0.1) distributed within HEASOFT. We 
used the latest release of the XRT Calibration Database 
(CALDB 2.6) and applied standard screening criteria to 
obtain the final cleaned event list. We selected only 
events with grades 0-12 for PC mode data and 0-2 for 
WT mode data. Such a selection provides the best com- 
bination of spectral resolution and detection efficiency. 
Our analysis has been performed over the 0.3-10 keV en- 
ergy band. 
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TABLE 1 

Results from X-ray lightcurve fits. Two different 

MODELS HAVE BEEN CONSIDERED AND THEY ARE DISCUSSED IN 
THE TEXT. 



Parameters 


Model 1 


Model 2 




Multiple broken PL 


Simple broken PL 



Fig. 2. — Upper panel: XRT light curve of GRB070110 in the 0.3-10 keV band. The two models described in the text are shown: power 
law segments with three temporal breaks (solid line) and a simple broken power law (dashed line). The bump at t~530 s is modeled with 
a Gaussian function, the late one at t~5xl0 4 s with a FRED profile. The four phases of the X-ray light curve are marked: (I) an early 
decay, (II) an apparent plateau followed by (III) a rapid drop, and (IV) a final shallow decay. Middle panel: Hardness ratio (H/S) light 
curve. It compares source counts in the hard band (H: 1-10 keV) and in the soft band (S: 0.3-1 keV). Lower panel: Photon index Txrt 
temporal variations. These values were derived fitting the X-ray spectra with an absorbed power law model. 

2.3.1. Temporal analysis 

We extracted WT data in a rectangular region, 40 x 20 
pixels wide, centered on the source position. The back- 
ground contribution was estimated from a region, with 
the same shape and size, sufficiently offset (>2 arcmin) 
from the source position to avoid contamination from the 
PSF wings and free from contamination by other sources. 

The first 300 s of PC observations were affected by pile- 
up. To account for this, we chose an annular extraction 
region centered on the source, having an inner radius of 4 
pixels and an outer radius of 30 pixels. From the second 
orbit on (t>3xl0 4 s), when pile-up is no longer present, 
the source count rate was estimated in a circular region 
with a 30 pixels radius. The count rate evolution in the 
later observations (t>3xf0 5 s) was obtained using the 
task sosta of the ximage package, using an extraction 
region that optimizes the signal to noise ratio of the de- 
tection. 

All the light curves presented here are background sub- 
tracted and corrected for Point Spread Function (PSF) 
losses, vignetting effects and exposure variations. 

The 0.3-10 keV light curve is shown in Fig. 2 (upper 
panel). WT and PC data were binned to achieve a min- 
imum signal to noise ratio of 8 and 5, respectively. Late 
points (t>10 5 s) are 3a detections. We also included the 
detection taken in Image (IM) mode, converting Digi- 
tal Number (DN) units into count rate (Hill et al. 2006; 
Mangano et al. 2007b). 

The light curve has been modeled with power law seg- 
ments of different slopes, whose best fit parameters are 
reported in the second column of Table 1. The best fit 
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model is shown by the solid line in the upper panel of 
Fig. 2. We also performed a fit of the light curve with 
a simple broken power law, excluding from the analysis 
the plateau and the following steep drop. The model is 
shown by the dashed line in Fig. 2, with best fit param- 
eters reported in Table 1 (third column). 

In both models we used a Gaussian function, centered 
at t~530 s and a ^50 s wide, to model the small bump at 
the end of the early decay. The fast rise, after the abrupt 
drop at ~2xl0 4 s, and the following decay are well de- 
scribed by a Fast-Rise-Exponential-Decay (FRED) pro- 
file, peaking at t~55 ks. 

The hardness ratio light curve is reported in the mid- 
dle panel of Fig. 2. It is defined as the ratio between the 
source counts in the hard band (H: 1-10 keV) and the 
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TABLE 2 

Results from X-ray spectral fits 3 -. 



X-ray 


Instr. 


Start 


Stop 


r x 


X 2 /dof 


Phase 


Mode 


(s) 


(<0 






Early 


W 1 


99.7 


11 1 


1 oa _l_n i 1 
l.o4±U.ll 


14/21 


decay 


W 1 


117 
11 I 


137 


1 no_l_n i a 

i.yz±u. 14 


23/17 




WT 


137 


177 


2.16±0.14 


21/18 




WT 


177 


265 


2.21±0.14 


6/14 




PC 


266 


644 


1.7±0.3 


7/6 


Plateau 


PC 


4043 


5237 


2.32±0.17 


12/15 




PC 


5237 


6426 


2.05±0.16 


21/17 




PC 


9825 


10997 


2.06±0.14 


23/16 




PC 


10997 


12210 


1.92±0.13 


13/16 




PC 


15607 


16896 


2.16±0.14 


28/19 




PC 


16896 


17993 


2.17±0.14 


20/17 


Drop 


PC 


21394 


29559 


2.12±0.15 


14/15 


Shallow 


PC 


5.8xlO b 


1.4xl0 b 


2.17±0.20 





dccay b 



a Spcctra were modelled with an absorbed power law. Galac- 
tic and intrinsic absorption were kept fixed at the value of 
1.86X10 20 cm~ 2 and 2.6X10 21 cm~ 2 , respectively. 

b The Cash statistic was applied. 



source counts in the soft band (S: 0.3-1 keV). An initial 
hard-to-soft evolution is present followed by a hardening 
of the spectrum, which corresponds to the small bump 
in the X-ray light curve at ~530 s. The hardness ra- 
tio shows an increasing trend at the beginning of the 
observed plateau and spectral variations throughout it. 
The very steep drop displays a hard-to-soft evolution. 
Later points (t>10 5 s) have a harder spectrum, maybe 
due to a late flaring activity, as the fluctuations in the 
light curve suggests. At times later than 2xl0 5 s the 
afterglow evolution does not show any further spectral 
variation. 

2.3.2. Spectral analysis 

In order to quantify the spectral variations seen in the 
hardness ratio light curve, we performed a time-resolved 
spectral analysis. 

Source and background spectra were extracted from 
the same regions used to create light curves (§ 2.3.1). 
The relevant ancillary response files were generated us- 
ing the task xrtmkarf . Time intervals were selected ac- 
cording to light curve phases and to have at least 400 
net counts each. Only the first PC orbit spectrum (from 
266 s to 644 s) has a lower statistical quality ( ^200 
source counts), due to the presence of pile- up. Spectral 
channels were grouped so to have at least 20 counts each. 
The x 2 statistic was applied. 

All the X-ray spectra can be modeled with an ab- 
sorbed power law. The Galactic absorption component 
was kept fixed at the value of 1.86xl0 20 cm -2 (Dickey 
& Lockman 1990); an additional rcdshifted absorption 
component, modeling the host intrinsic absorption, was 
also included. In order to estimate the host absorp- 
tion we extracted a WT (from 100 s to 265 s) and a 
PC spectrum (from 4 ks to 30 ks). Since at low ener- 
gies (~0.5 keV) XRT spectra may be affected by calibra- 
tion uncertainties 1 , we performed our analysis excluding 
energy channels in the range 0.45-0.55 keV. From the 
joint spectral fit we obtained an intrinsic absorption of 

1 http://heasarc.gsfc.nasa.gov/docs/swift/analysis/ 



TABLE 3 

Swift UVOT photometry of GRB 070110 
afterglow. 



Start 


At 


Magnitude 3 


Error b 


Filter 


(s) 


(s) 








103 


100 


19.90 


0.31 


White 


4249 


200 


19.18 


0.26 


White 


5681 


200 


19.69 


0.18 


White 


85 


9 


17.35 


UL 


V 


208 


400 


20.28 


0.45 


V 


4658 


200 


18.94 


0.33 


V 


6091 


200 


18.83 


UL 


V 


22298 


18665 


20.38 


0.27 


V 


161578 


47181 


21.52 


0.32 


V 


4044 


1632 


20.33 


0.30 


B 


11644 


598 


20.59 


0.21 


B 


28983 


603 


21.32 


0.36 


B 


114817 


12775 


22.08 


0.32 


B 


662 


15 


19.23 


0.51 


U 


5272 


200 


19.63 


0.22 


u 


10733 


906 


19.84 


0.14 


u 


17423 


602 


20.22 


0.23 


u 


28072 


906 


20.42 


0.22 


u 


34762 


603 


20.84 


0.31 


u 


235449 


135190 


22.58 


UL C 


u 


495668 


81911 


22.07 


UL 


u 



a Values are not corrected for extinction. 
Errors at the 68% confidence level or 3(7 upper limits 
(ULs) arc given. 

c Our re-analysis using the refined optical position 
gives a 3(7 UL instead of the detection quoted in 
Krimm et al. (2007a). 

N h f st = (2.6±l.l)xl0 21 cm" 2 . This value was derived as- 
suming a host galaxy with solar metallicity. 

We then performed time-resolved spectral fits keeping 
the absorption fixed as above and leaving only the photon 
index Tx and the normalization as free parameters. We 
also tested whether the initial softening of the spectrum 
could be due to a decreasing intrinsic absorption, as pre- 
viously witnessed in other GRBs (Cusumano et al. 2007, 
Campana et al. 2006), but the low number of counts did 
not allow us to constrain the Ng ost behavior. 

In order to estimate the spectral index during the fi- 
nal shallow decay, we extracted a spectrum in the time 
interval from 5.8xl0 5 s to 1.4xl0 6 s, using a circular re- 
gion of 5 pixels radius centered on the source position. 
Because of the low number of counts (~120) and the neg- 
ligible background contamination, the Cash statistic was 
applied (Cash 1979). We obtained r x =2.17±0.20. 

The selected time intervals and results from the spec- 
tral fit arc listed in Table 2. The photon index variations 
with time are reported in the lower panel of Fig. 2. 

2.4. Ultraviolet/ Optical data 

Independent measurements of the position of the af- 
terglow were made after summing the images taken dur- 
ing the first day after the trigger for the White, V, B, 
and U filters. The mean position is R.A.=00 h 03 m 39.23 s , 
Dec.=-52°58'26.9" with an estimated lcr uncertainty of 
0.25" in each direction, based on the rms deviation of the 
four measurements. This position is consistent with the 
ground-based position given in Malesani et al. (2007b) . 

We performed photometry using a circular aperture of 
radius 2.5" centered on the position of the optical after- 
glow. This choice allows us to minimize the background 



Swift observations of GRB070110 



5 



o 
o 
o 



o 
o 



-3 O 

=1 



c 

CD 
T3 

X 
13 



O 

d 



i i 1 1 1 in 1 — i i 1 1 1 1 n — ri — i i 1 1 1 1 1| 1 — i i 1 1 1 1 1 1 1 — i i i.i i in — n — i i 1 1 1 ii 



1 1 I I I I I l| 1 TTT] 

+ BAT @ 1 keV 

+ XRT @ 1 keV 

• U filter 

* B filter 
■ V filter 
x White filter 




10 



100 



1000 



10 4 



10 5 



10 6 



Time since BAT trigger (s) 



Fig. 3. — Combined BAT, XRT and UVOT light curves. BAT and XRT count rates were converted into monochromatic fluxes using the 
results from the spectral analysis. UVOT magnitudes were corrected for the Galactic extinction and converted into flux densities at the 
central wavelength of each filter. In the plot UVOT flux values have been scaled by a factor 10 — 3 . The optical best fit models, discussed 
in the text, are shown by solid lines. The dashed line shows the broken power law model (model 2). The vertical dashed-dotted lines mark 
the times at which SEDs were computed. 



contribution. Aperture corrections were computed us- 
ing stars in the images to convert the 2.5" photometry 
to a 6" aperture for the V, B and U filters, and a f2" 
aperture for the White filter. Adjacent exposures with 
AT/T<0.1 were coadded to improve the signal to noise. 

The instrumental magnitudes were transformed to 
Vega magnitudes using the standard photometric zero 
points in the Swift /UVOT calibration database 2 . Ta- 
ble 3 contains the results of the UVOT photometry. The 
reported errors and upper limits (ULs) are at ler and 3cr 
confidence level, respectively. 

The Galactic reddening in the burst direction is E(B- 
V) = 0.014 (Schlegel et al. 1998). Using the value from 
the extinction curve in Pei (1992) for the Milky Way at 
the central wavelength of each filter, we estimated the 
extinction for the UVOT filters to be Ay = 0.04 mag, 
Ab = 0.06 mag and Au = 0.07 mag. 

Fig. 3 presents the multi-wavelength light curve of 
GRB 070110, including 7-ray, X-ray and UV/optical 
data. Both BAT and XRT light curves were converted 
into flux units using the best fit spectral models. The 
V, B and U magnitudes were corrected for Galactic ex- 
tinction along the line of sight and then converted to 
monochromatic fluxes at the central wavelength of each 
filter. Since the afterglow is not detected blueward of the 
U filter, we took 4450A as the effective mid-wavelength 

2 http: / /heasarc.gsfc. nasa. gov / docs /heasarc /caldb /swift /docs/ uvot / i: 



of the White filter. The conversion factor from count rate 
to monochromatic flux density depends on the afterglow 
spectral shape and the White filter effective area curve 
and it was found to be 0.026 mJy s cts -1 . 

With the exception of the first detection in the U, 
V and White filter, the UVOT light curves show a de- 
caying behavior that can be well described by a simple 
power law. An independent fit of the three light curves 
gives the following temporal slopes: ay=0.65±0.14, 
a B =0.53±0.13, at7=0.54±0.15. We note that the UVOT 
late slopes are consistent, within the errors, to the late- 
time X-ray slope in Table 1. Including in the fit also 
the late VLT detection in the V band, reported in Male- 
sani et al. (2007a), the slope ay steepens to 0.88±0.09. 
Alternatively, we attempted a fit with a broken power 
law, which gives an initial decay of 0.62±0.16, steepen- 
ing at ~150 ks to 1.3±0.2. Both models provide a good 
description of the light curve in the V band. 

A joint fit of the three UVOT light curves (U, B, and 
V) has been performed using the functional form of Will- 
ingale et al. (2007). We allowed the normalizations to 
vary independently, tying other parameters. Upper lim- 
its and the VLT detection were not considered in the fit. 
We found that our light curves follow a power law decay 
with a common slope of 0.63±0.13. 

The best fit models of the three UVOT light curves 
(U, B, and V) with the Willingale et al. (2007) function 
kheWn by solid lines in Fig. 3. The fast rise of the 
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optical emission is poorly sampled and we note that the 
first detection in the White filter suggest a more gentle 
onset of the afterglow emission. 

For a comparison we report the broken power law 
model (model 2), described in § 2.3.1 (dashed line in 
Fig. 3). 

2.5. Optical/X-ray Spectral Energy Distributions 

We calculated the Spectral Energy Distribution (SED) 
at four representative times: t=400 s during the initial 
decay, t=5.5 ks and t=17 ks at the beginning and at 
the end of the plateau phase, and t=100 ks during the 
shallower decay. The four times are marked by vertical 
lines in Fig. 3. 

Fig. 4 shows the four SEDs derived from the optical 
and X-ray light curves. Flux values in the UVOT filters 
at the selected times are corrected for Galactic extinc- 
tion (empty symbols) and were derived from the best fit 
models shown in Fig. 3. At a redshift of z=2.352 the 
Lya forest is present throughout the U filter. Therefore, 
in constructing the SEDs we have corrected the U band 
count rate for the expected intergalactic Lya transmis- 
sion of 86% over the U passband (Madau 1995). 

X-ray data were converted into flux units for each time 
selecting the appropriate photon index value from Ta- 
ble 2. X-ray fluxes and their error regions are shown by 
the black cones. The contribution of the host galaxy to 
extinction have been estimated as the additional extinc- 
tion required to have optical and X-ray data lying on the 
same power law at t=100 ks. Accounting for a host ex- 
tinction of Av=0-0S mag with the extinction curve for 
the Small Magellanic Cloud (SMC; Pci 1992) provides a 
better match between optical and X-ray data. The entire 
data set at the later times seems consistent with a com- 
mon physical origin. UVOT fluxes corrected for both the 
Galactic and the intrinsic extinctions are shown by filled 
symbols. The best fit power law of the optical SED is 
shown by the dot-dashed line. 

At early times (t=400 s) the optical data are not con- 
sistent with the extrapolation of the X-ray spectrum to 
low energies. At 5.5 ks and 17 ks, during the apparent 
plateau, the optical and X-ray spectral distributions are 
also completely inconsistent with one another, implying 
different origins for the optical and X-ray photons. 

In § 2.3.1 we have shown that the X-ray light curve, 
excluding the plateau, can be fitted with a simple two 
component broken power law plus flares. Since the slope 
of this power law is marginally consistent with that ob- 
served in the optical bands, it follows that the SED con- 
sistency observed at late times also applies to this com- 
ponent. This is illustrated in Fig. 4 (top right panel and 
bottom left panel) , where we have extrapolated the late- 
time SED to 5.5 and 17 ks, using light curve model 2 
(third column in Table 1) to evaluate the X-ray count 
rates. These values were derived using the photon index 
value at late times (t>580 ks). The extrapolated values 
and their error regions are shown by the grey cones. 

If we assume for the host galaxy of GRB 070110 
the same gas-to-dust ratio measured in the SMC, 
Nh/^4v=1-6x10 22 cm~ 2 mag -1 (Wcingartner & Draine 
2000), an intrinsic hydrogen column density of 
1.3xl0 21 cm~ 2 is inferred from the estimated rest frame 
visual extinction Ay=0.08 mag. 

The intrinsic absorption derived from fitting the X-ray 



spectra is instead (1.6±0.7) x 10 22 cm" 2 for a SMC-like 
absorber, higher than the one expected from the gas- 
to-dust relationship. We derived the latter absorption 
value modeling the intrinsic X-ray absorption spectral 
component with the model tbvarabs in XSPEC (Arnaud 
1996), keeping the absorber abundances fixed to SMC 
values (~Z©/8). 

The high X-ray absorption column and low optical ex- 
tinction measured for GRB 070110 require a higher than 
expected gas-to-dust ratio, and are indicative of an in- 
terstellar medium significantly different from the one ob- 
served in the local universe. Such difference could be 
ascribed to the action of the GRB radiation on its imme- 
diate environment (Waxman & Draine 2000; Pcrna et al. 
2003) or, alternatively, to an intrinsically higher gas- 
to-dust ratio in GRBs birth sites (Watson et al. 2006). 
Discrepancies between optical and X-ray extinctions in 
GRBs environments were first found by Galama & Wijers 
(2001) and then confirmed by following studies on large 
GRBs samples (e.g. Kann et al. 2006; Starling et al. 
2006; Tagliaferri et al. 2007). 

3. DISCUSSION 
3.1. The Early X-Ray Tail 

As shown in Fig. 2 (upper panel), the X-ray light curve 
is composed of four distinct components. It starts with 
(I) an early steep decay component with a = 2.44 ±0.13 
following the prompt gamma-rays. The early decay is 
followed by (II) a plateau with an essentially constant 
flux which extends up to t^20 ks. Luminosity fluctua- 
tions are seen throughout the plateau phase. Following 
the plateau is (III) a remarkable steep drop with a tem- 
poral index of ~9.0. After this steep fall the afterglow 
light curve rises again, showing a late flare, then enters 
(IV) a more normal decay segment with a ~0.7. 

The hardness ratio light curve (Fig. 2, middle panel) 
shows an initial hard-to-soft spectral evolution during the 
early decay. Our spectral analysis confirms a softening 
of the spectral index /3x from a value of 0.8 to 1.2. The 
early X-ray decay has been generally interpreted as the 
prompt emission tail, due to the delay of propagation of 
photons from high latitudes with respect to the line of 
sight (Tagliaferri et al. 2005; Nousek et al. 2006; O'Brien 
et al. 2006; Zhang et al. 2006a). 

The curvature effect alone, which suggests a=2+/3 
(Kumar & Panaitescu 2000), cannot explain both the 
spectral evolution feature and the shallower than ex- 
pected temporal slope. Such behavior has been seen in 
other GRBs detected by Swift and, as reported by Zhang 
et al. (2006b) , spectral evolution seems a common feature 
in early, bright GRB tails (e.g. GRB 060614, Mangano 
et al. 2007a). Phenomenological models involving an un- 
derlying central engine afterglow with a steep spectral 
index or an internal shock afterglow with a cooling fre- 
quency passing through the XRT window may be candi- 
dates to interpret these spectral evolutions (Zhang et al. 
2006b). 

3.2. The Plateau Phase 

The observed plateau is a feature of great interest. 
This component displays an apparently constant inten- 
sity extending up to t^20 ks, followed by an abrupt 
drop with a very steep decay index. Compared with 
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Fig. 4. — Optical, ultraviolet and X-ray SEDs at 400 s, 5.5 ks, 17 ks, and 100 ks. UVOT points have been corrected for the Galactic 
extinction Av=0.04 mag (empty symbols) and for the host extinction Av=0.08 mag (filled symbols) along the line of sight . The dashed 
line is the best fit of the optical data corrected only for Galactic extinction. The dot-dashed line is the best fit of the optical data corrected 
for both the Galactic and the host intrinsic extinctions. The solid black lines define the cone corresponding to the 90% uncertainty on the 
spectral slope in the X-ray band. This cone is extrapolated to the optical band through the dotted lines. The X-ray fluxes, extrapolated 
at t=5.5 ks and t=17 ks according to our model 2. and the corresponding error cones are also shown (grey solid lines). 



the canonical XRT light curves GRBs observed by Swift 
(Nousek et al. 2006; O'Brien et al. 2006; Zhang et al. 
2006a; Willingale et al. 2007), such a steep decay follow- 
ing a plateau is unique. 

In most other Swift GRBs, the plateau is followed by 
a "normal" decay which could be generally interpreted 
as the standard external forward shock afterglow. The 
plateau in those cases is therefore consistent with a re- 
freshed shock (Rees & Meszaros 1998; Zhang et al. 2006a; 
Nousek et al. 2006) as the total energy in the blast-wave 
is replenished with time. 

The abrupt steep decay in GRB 070110, on the other 
hand, is inconsistent with an external shock origin of 
the plateau. This is because, within the external shock 
model, it is impossible to produce a decay slope as steep 



that observed. The steepest possible decay from an ex- 
ternal shock occurs after a jet break. The post jet break 
decay slope can be as steep as a ~ p, where p is the elec- 
tron spectral index (Rhoads 1999; Sari et al. 1999) and 
afterglow modeling suggests that p is typically between 
1.5-3. 

A steep decay may be expected when the blast-wave 
leaves a density clump and falls into a very tenuous 
medium. In such a case, the curvature effect (Kumar & 
Panaitescu 2000) defines the decay slope. However, de- 
tailed numerical calculations suggest that the zero-time 
point of the external shock scenarios cannot substantially 
be shifted from the GRB trigger time (Zhang et al. 2006a; 
Lazzati & Begelman 2006; Kobayashi & Zhang 2007). 
As a result, the expected steep decay index within this 
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Fig. 5.— Combined BAT and XRT light curves of GRB 050904 
(2=6.29) and GRB 070110 (2=2.35) in the source rest frame. 

scenario is a~2-r-/?~3, far shallower than the observed 
value. 

By assuming that the steep decay following the plateau 
is due to the curvature effect, we performed the same test 
of Liang et al. (2006), who tested the internal origin of 
X-ray flares. We found that to satisfy the a — 2 + j3 
condition the zero time to needs to lie in the range 16- 
18 ks after the trigger, slightly before the beginning of the 
abrupt drop. The very steep decay component therefore 
strongly suggests that the plateau is of internal origin, 
since to is allowed to shift to much later epochs within 
the internal scenarios (Zhang et al. 2006a; Liang et al. 
2006). 

Several other GRBs have previously shown similar, 
but not exactly the same, behavior. The afterglow of 
GRB 051117A (Goad et al. 2006) shows a sharp drop 
of the X-ray flux at ~10 4 s, followed by a shallower de- 
cay (a~0.7). However, in this case the early X-ray light 
curve seems dominated by flaring activity and it does 
not show any evident plateau. Fig. 5 compares the rest 
frame light curves of both GRB 050904 and GRB 070110. 
GRB 050904 (Cusumano et al. 2006, 2007) shows a re- 
markably similar behavior to GRB 070110, with a rapid 
decay following the multiple episodic X-ray flares. How- 
ever, GRB 070110 displays a much smoother emission 
episode before the rapid decay. While the steep decay fol- 
lowing the flares in GRB 050904 is expected, the dropoff 
after the plateau in GRB 070110 is not. 

On the basis of Swift data, collected in portions last- 
ing less than an hour, we cannot definitely rule out 
that the observed plateau is due to the same mechanism 
producing erratic X-ray flares in other GRBs, such as 
GRB 050904. The hardness ratio (Fig. 2, middle panel) 
appears to rise slowly at the start of the plateau (~4 ks) 
and a decreasing trend is visible during the sharp de- 
cay (~20 ks). Such behavior is similar to the one seen 
in the hardness ratio light curves of X-ray flares, al- 
though timescales associated with a single flare are usu- 
ally shorter. Attributing the whole emission between 4 ks 
and 20 ks to a superposition of many flares allows the 
final drop to be explained with the curvature effect, re- 
setting the reference time to just before the rapid de- 
cay. In this case, the smoothness of the flat phase could 
be a selection effect, due to lacking full coverage of the 
light curve. The hypothesis of flares blending to form an 



apparent plateau seems however very unlikely, requiring 
a high degree of homogeneity between successive X-ray 
flares, since no large amplitude variations have been ob- 
served (see Fig. 2). 

We conclude that GRB 070110 robustly suggests new 
properties of the GRB central engine: it must be able 
to last for an extended period of time (in this case up 
to 8xl0 3 s in the source rest frame) with an essentially 
constant radiation power. 

Previous Swift observations of X-ray flares (Burrows 
et al. 2005b; Barthelmy et al. 2005b; Falcone et al. 2006; 
Romano et al. 2006; Campana et al. 2006; Cusumano 
et al. 2006, 2007) have revealed that the GRB central en- 
gine can indeed be active to such a late epoch. However, 
the energy output (at least the observed luminosity) in 
previous cases has to be intermittent. The mechanisms 
to interpret these intermittent X-ray flares therefore re- 
quire that the central engine accretes episodically, either 
due to fragmentation or density fluctuation of the accre- 
tion disc (e.g. King et al. 2005; Perna et al. 2006; Ross- 
wog 2006), or due to modulation of the magnetic fields 
near the accretor (Proga & Zhang 2006). Alternatively, 
the intermittent behavior of the central engine may be 
related to the unsteady magnetic activity of the central 
engine (e.g. Dai et al. 2006). None of these sugges- 
tions seems to be able to interpret the internal plateau in 
GRB 070110 straightforwardly. Rather we require that 
the central engine is continuously active for a long period 
of time. 

Such long-term central engine activity has been dis- 
cussed before (e.g. Zhang & Meszaros 2001 for a gen- 
eral discussion of a central engine with decaying lumi- 
nosity L(t)oct~ g ). The most straightforward example is a 
spinning-down pulsar (Dai & Lu 1998; Zhang & Meszaros 
2001). This model has been introduced to interpret the 
external plateau as due to refreshed shocks, although 
those data can be also interpreted without introducing 
such a long-term central engine (Rees & Meszaros 1998; 
Granot & Kumar 2006). The discovery of the internal 
shock origin of the plateau in GRB 070110 breaks the de- 
generacy of the refreshed shock interpretation, and gives 
the first direct evidence of a continuous long-lasting cen- 
tral engine. 

We propose that the engine powering the plateau could 
be a spinning-down pulsar, which has a constant lumi- 
nosity lasting for an extended period of time (Shapiro 
& Teukolsky 1983; Zhang & Meszaros 2001). The du- 
ration of the plateau depends on the unknown pulsar 
parameters, but given a reasonable radiation efficiency, 
the luminosity (~10 48 ergs -1 ) and the observed duration 
(~16 ks) of the plateau are consistent with the parame- 
ters of a new-born magnetized millisecond pulsar as the 
GRB central engine. 

According to Zhang & Meszaros (2001), the continu- 
ous injection luminosity L C m.o and the charateristic time 
scale T em , when the plateau breaks down, are related to 
the pulsar initial parameters: 

L em>0 ~ 10 49 B 2 p , 15 PoT-3^6 erg s" 1 . (1) 

T om ~ 2.05 x 10 3 h 5 B- 2 15 P^ 3 Rs 6 s. (2) 

where B p = B p i 5 x 10 15 G is the dipolar field strength at 
the poles, Po,-3 is the initial rotation period in millisec- 
onds, J45 is the moment of inertia in units of 10 45 gem 2 , 
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and i?6 is stellar radius in units of 10 6 cm. 

If we assume that a significant fraction of the spin- 
down luminosity is emitted in the X-ray band and take 
standard values of I i5 ~ 1 and i? 6 ~ 1, from the 
plateau parameters we can infer a pulsar initial period 
of P a <1 ms and a magnetic field B p >3x 10 14 G. The 
energy of the plateau puts only a lower limit to the real 
energy of the central pulsar L cnlj o and without a good 
estimate of the radiation efficiency it is not possible to 
constrain better the pulsar parameters. It is worth to 
note that our estimate of Po is very close to the break- 
up period for a neutron star, P m i n =0.96 ms (Lattimcr 
& Prakash 2004), so the derived pulsar parameters can 
be taken as good approximations of the right values. 

The observed decay slope following the plateau is much 
steeper than the model prediction (a<~2). However, ac- 
counting for detailed energy dissipation mechanisms and 
possible magnetic field decay at the central engine can 
steepen the decay (Zhang et al., in preparation). The 
energy dissipation mechanism of a spinning-down pulsar 
before deceleration is not well studied, but it may be 
related to the breakdown of the magneto-hydrodynamic 
condition in the magnetized outflow (Usov 1994; Spruit 
et al. 2001; Zhang & Meszaros 2002a) and powered by 
magnetic reconnections (Drcnkhahn & Spruit 2002). In- 
stability inside the energy dissipation regions would in- 
duce fluctuations in the emission region, giving rise to the 
flickering light curve and the hardness ratio variations on 
the plateau. 

3.3. The Final Shallow Decay 

The late X-ray light curve shows a shallow power law 
decay, on which a flare at ^55 ks and mini-flares up 
to ^200 ks are superimposed. These flares could be 
again interpreted as late central engine activities (Bur- 
rows et al. 2005b; Romano et al. 2006; Falcone et al. 2006; 
Zhang et al. 2006a; Fan & Wei 2005) which are likely 
related to the episodic accretion processes and may be 
different from the spindown-powered plateau. 

From the SEDs analysis (§ 2.5) we derived that in the 
time interval 5-20 ks we observe an internal afterglow 
component, which contributes strongly to the X-ray but 
not the optical band, and a second afterglow compo- 
nent, likely of external shock origin, which tracks the 
optical light curves. This component dominates after 
the late steep drop, when the emission from the central 
engine switched off. The late optical/X-ray spectrum 
(t=100 ks) can be described by a continuous power law 
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(/3=1.00±0.14), indicating that the optical and the late 
X-ray afterglow may arise from the same physical com- 
ponent. 

The shallow temporal slopes in the X-ray band 
(ax ~0.7) and in the optical bands (a op t=0.63±0.13) are 
consistent. However we note that the observed temporal 
decays are not in agreement with the standard closure 
relations (e.g. Zhang & Meszaros 2004), which predict 
a temporal slope a >1 for both the X-ray and the opti- 
cal light curves. Such shallow behavior at late times is 
seen in ~50% of the X-ray light curves observed by Swift, 
as noted by Willingale et al. (2007), and it may suggest 
that the external shock is kept refreshed by a long-lasting 
energy injection. 

4. CONCLUSION 

GRB 070110 is a long burst (T 90 ~90 s) with a standard 
behavior in the 7-ray band, that shows extraordinary 
properties of the X-ray afterglow. The presence of a 
plateau followed by an unexpected steep decay implies 
a new emission component from a long-lasting central 
engine. Such constant emission may be generated by a 
spinning down pulsar. This mechanism indeed allows us 
a better interpretation of the GRB070110 phenomenol- 
ogy than the superposition of multiple flares. However 
the naive picture we discussed cannot account for the 
rapid decay at the end of the plateau, indicating that a 
more detailed theoretical study is needed. A drop of the 
radiation efficiency, caused by the pulsar deceleration, 
or a decaying magnetic field of the new-born system 
could help to explain the steep decay rate. 
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